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Abstract 

We estimate the two-photon exchange contribution to elastic electron-proton scattering at large 
momentum transfer by using a quark-parton representation of virtual Compton scattering. We 
thus can relate the two-photon exchange amplitude to the generalized parton distributions which 
also enter in other wide angle scattering processes. We find that the interference of one- and 
two-photon exchange contribution is able to substantially resolve the difference between electric 
form factor measurements from Rosenbluth and polarization transfer experiments. Two-photon 
exchange has additional consequences which could be experimentally observed, including nonzero 
polarization effects and a positron-proton/electron-proton scattering asymmetry. The predicted 
Rosenbluth plot is no longer precisely linear; it acquires a measurable curvature, particularly at 
large laboratory angle. 
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I. INTRODUCTION 



There are two experimental methods for extracting the ratio of electric (G P E ) to mag- 
netic (G P M ) proton form factors from electron-proton scattering: unpolarized measurements 
employing the Rosenbluth separation technique, and polarization experiments. In the lat- 
ter case, one measures the correlation of the spin of the incident polarized electron with 
the polarization components of the outgoing proton, parallel Pi or perpendicular P s (in the 
scattering plane) to its momentum . The ratio of cross sections for the two outgoing 

proton polarizations gives Ge/Gm directly: 



P l \jT(l + e) G M (Q 2 ) 
The kinematic functions e and r are 



P. _ I 2e G E (Q 2 ) (1) 



and 



- = 1 + 2(1 + r)tan 2 ^ , (3) 

£ 2 

where Q 2 = —q 2 = —t is the momentum transfer squared, 9 is the laboratory scattering 
angle, and < e < 1. Equivalent information may be obtained in scattering of longitudinally 
polarized electrons on a polarized proton target. 

The Rosenbluth method relies on measuring the differential cross section 

« G 2 M + E -G\ , (4) 



dflLab T 

with the proportionality factor being well known, and isolating the e dependent term. In 
each case, the extraction method for Ge/Gm assumes single-photon exchange between the 
electron and nucleon. 

Recent polarization experiments at the Thomas Jefferson Laboratory (JLab) have 
confirmed the earlier Rosenbluth measurements from SLAC [6]. However, at large Q 2 , all of 
the Rosenbluth measurements are at distinct variance with JLab measurements of G p ^/G p , r 

nn 

obtained using the polarization technique |l|, |2| . Since G P E contributes to the unpolarized 
cross section at only a few percent level for the Q 2 range in question, it is necessary to 



2 



identify any possible systematic corrections to the Rosenbluth measurements at the percent 
level which could be responsible for this discrepancy. 

One possible explanation for the discrepancy between the Rosenbluth and polarization 
methods is the presence of two-photon exchange effects, beyond those which have already 
been accounted for in the standard treatment of radiative corrections. A general study 
of two- (and multi)-photon exchange contributions to the elastic electron-proton scattering 
observables was given in In that work, it was noted that the interference of the two- 
photon exchange amplitude with the one-photon exchange amplitude could be comparable 
in size to the {G P E ) 2 term in the unpolarized cross section at large Q 2 . In contrast, the 
two-photon exchange effects do not impact the polarization-transfer extraction of Ge/Gm 
in an equally significant way. Thus a missing and unfactorizable part of the two-photon 
exchange amplitude at the level of a few percent may well explain the discrepancy between 
the two methods. 

Realistic calculations of elastic electron-nucleon scattering beyond the Born approxima- 
tion are required in order to demonstrate in a quantitative way that two-photon exchange 
effects are indeed able to resolve this discrepancy. In particular, one wants to study quan- 
titatively the "hard" corrections which will arise when both exchanged photons are far off 
shell or the intermediate nucleon state suffers inelastic excitations. Calculations of these 
corrections require a knowledge of the internal structure of the nucleon and thus could not 
be induded in the da.ie Eff and were not Included in the m0r e recent HQ e al cn la t 10 n S 
of radiative corrections to eN elastic scatte ring . 

A first step was performed recently in [12|, where the contribution to the two-photon 
exchange amplitude was calculated for the elastic nucleon intermediate state. In that calcu- 
lation it was found that the two-photon exchange correction with an intermediate nucleon 
has the proper sign and magnitude to partially resolve the discrepancy between the two 
experimental techniques. 

In an earlier short note Q], we reported the first calculation of the hard two-photon 
elastic electron-nucleon scattering amplitude at large momentum transfers by relating the 
required virtual Compton process on the nucleon to generalized parton distributions (GPD's) 
which also enter in other wide angle scattering processes. This approach effectively sums all 
possible excitations of inelastic nucleon intermediate states. We found that the two-photon 
corrections to the Rosenbluth process indeed can substantially reconcile the two ways of 
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measuring Ge/Gm- Our goal in this paper is to give a detailed account of our work, and to 
present numerical results for a number of quantities not included in the shorter report. 

Perturbative QCD factorization methods for hard exclusive processes provide a systematic 
method for computing the scaling and angular dependence of real and virtual Compton 
scattering at large t. For example, PQCD predicts that the leading-twist amplitude for 
Compton scattering 7^ — > ^yp can be factorized as a product of hard-scattering amplitudes 
TnilQQQ — * where the quarks in each proton are collinear, convoluted with the initial 

and final proton distribution amplitudes <f)(xi,Q) Q|. All of the hard-scattering diagrams 
fall at the same rate at large momentum transfer whether or not the photons interact on 
the same line. Although, the predictions for the power-law falloff and angular dependence 
of Compton scattering are consistent with experiment, the leading-twist PQCD calculations 
of the wide angle Compton amplitudes appear to substantially underpredict the magnitude 
of the observed Compton cross sections |l5| . 

Since an exact QCD analysis of virtual Compton scattering does not appear practical, 
we have modeled the hard two-photon exchange amplitude using the "handbag approxi- 
mation" Q], in which both photons interact with the same quark. The struck quark is 
treated as quasi-on shell. In particular, we have neglected the amplitudes where the two 
hard protons connect to different quarks, the "cat's ears" diagrams, as well as the diagrams 
in which gluons interact on the fermion line between the two currents. The handbag di- 
agrams contain the "J = 0" fixed pole, the essential energy-independent contribution to 
the real part of Compton amplitude which arises due to the local structure of the quark 
current The handbag approximation has proven phenomenologically successful in 

describing wide-angle Compton scattering at moderate energies and momentum transfers. 
As we shall show, the handbag approximation allows the two-photon exchange amplitude 
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y 



2l| . thus providing 



to be linked to the generalized parton distributions (GPD's) 
considerable phenomeno logical guidance. 

Brooks and Dixon and Vanderhaeghen et al. have shown that PQCD diagrams where 
the photons attach to the same quark dominate the Compton amplitude on the proton, 
except at backward center-of-mass angles. 1 . The dominance of the handbag diagrams in the 



1 This suppression of the cat's ears diagrams at forward angles could be due to the momentum mismatches 
which occur when photons couple to different quarks. Another possible explanation is that in some 
kinematic regions, the cat's ears and handbag amplitudes have the same magnitude (or nearly so) except 
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the handbag approximation, 
have shown that electron- 



PQCD analysis provides some justification for the the use o 
However, it should be noted that Gunion and Blankenbecler 
deuteron scattering is dominated by the cat's ears diagrams at large momentum transfer 
provided that the deuteron wave function has Gaussian fall-off. The dominance of the 
handbag diagrams thus depends on the nature of the QCD wave functions, and the precise 
situation in the present case remains a subject for future study. 

Recently, a new category of Rosenbluth data has become available where the recoiling 
proton is detected The new data appear to confirm the older data, where the scattered 
electron was detected. The two-photon exchange contributions are the same whatever par- 
ticle is detected. However, the bremsstrahlung corrections, which are added to obtain an 
infrared finite result, are different. We shall defer detailed discussion of the proton-detected 
data until we can reevaluate the original proton-observed electron-proton bremsstrahlung 

2m as well as examine the radiative corrections which have been 



interference calculations 
applied to the new data 



123, l2£ 



The plan of this paper is as follows: 

The next section is devoted to kinematics, including the definitions of the invariants which 
define the scattering amplitudes and the formulae for the cross sections and polarizations 
in terms of those invariants. There are choices in the definitions of the invariants. We have 
presented the bulk of the paper with one choice; a sometimes useful alternative choice is 
summarized with cross section and polarization formulas in Appendix El Section ITTTI gives 
analytic results for the two-photon exchange scattering amplitudes at the electron-quark 
level, the hard scattering amplitudes required for the partonic calculation of two-photon 
exchange in electron-nucleon scattering. We have generally treated the quarks as massless. 
A quantitative discussion of modifications following from finite quark mass appears in Ap- 
pendix El Section ITVl details the embedding of the partonic amplitude within the nucleon 
scattering amplitude, using dominance of handbag amplitudes and GPD's. This section also 



for the charge factors. In these regions, the Compton amplitude would be proportional to the total charge 
squared (2e u + e^) 2 of the target proton. This is precisely the case in the low energy limit, where the 
Compton amplitude is indeed proportional to (2e u + ed) 2 = 1 for a proton. The result is reproduced by the 
handbag diagrams alone since, coincidcntally, 2e 2 + e 2 , = 1. In this scenario, the handbag approximation 
will fail for Compton scattering on a neutron or deuteron target. At higher energies, discussion of this 
scenario pertains to large angle Compton scattering, since in the forward direction the handbag diagrams 
are known to dominate. 
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discusses the particular GPD's which we have used in our numerical calculations. Section IVl 
shows numerical results, given graphically, for cross sections, single spin asymmetries, po- 
larization transfers, and positron-proton vs. electron-proton comparisons. Section also 
includes commentary about the possibility of extending the calculations to backward scat- 
tering (small values of |w|), and an assessment of how well two-photon physics reconciles the 
Rosenbluth and polarization transfer measurements of Ge/Gm- Section fVTl summarizes our 
conclusions. 



II. ELASTIC ELECTRON-NUCLEON SCATTERING OBSERVABLES 

In order to describe elastic electron-nucleon scattering, 

l(k, h) + N(p, X N ) - l(k', h') + N(p', X' N ), (5) 

where h, h', Xn, and X' N are helicities, we adopt the definitions 

p + p' k + k' , 

P = , K = — ^— , q = k-k =p -p , (6) 

define the Mandelstam variables 

s = {p + k)\ t = q 2 = -Q\ u = {p-k')\ (7) 

let v = K • P, and let M be the nucleon mass. 
The T-matrix helicity amplitudes are given by 

T^\ n = (k',h';p\X' N \T\k,h;p,X N ) . (8) 

Parity invariance reduces the number of independent helicity amplitudes from 16 to 8. Time 
reversal invariance further reduces the number to 6 Further still, in a gauge theory 

lepton helicity is conserved to all orders in perturbation theory when the lepton mass is zero. 
We shall neglect the lepton mass. This finally reduces the number of independent helicity 
amplitudes to 3, which one may for example choose as 

T +,+ . T _ +; +. T _ +; + = T +,+ _ (g) 

(The phase in the last equality is for particle momenta in the xz plane, and is valid whether 
we are in the center-of-mass frame, the Breit frame, or the symmetric frame to be defined 
below.) 
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Alternatively, one can expand in terms of a set of three independent Lorentz structures, 
multiplied by three generalized form factors. Only vector or axial vector lepton currents 
can appear in order to ensure lepton helicity conservation. A possible T-matrix expansion 
is (removing an overall energy-momentum conserving 5-function), 



p 2 f r 

^ G > _ t— F > 
M 



T h ,\'x N = 7^ < u(k', /i)7^u(fc, h) x u(p , X N ) 



u(p,X 



N 



u(k\ h)^ 5 u(k : h) x rz(p, \' N ) Yl 5 G' A u(p, X N ) \ . (10) 



Q 2 



This expansion is general. The overall factors and the notations G' M and F' 2 have been 
chosen to have a straightforward connection to the standard form factors in the one-photon 
exchange limit. 

There is no lowest order axial vector vertex in QED: the effective axial vertex in the 
expansion arises from multiple photon exchanges and vanishes in the one-photon exchange 
limit. One may eliminate the axial-like term using the identity 

s — u 

u(k')~/ ■ Pu(k) x u(p')^f ■ Ku(p) = — - — u{k')^^u(k) x u(p')^u(p) 

+ -u(k')^j 5 u(k) x u(p')Yl 5 u(p) , (11) 

which is valid for massless leptons and any nucleon mass. Hence, an equivalent T-matrix 
expansion is 

e 2 _ 

T h ,y N \ N = Q2 U ( k ,h)^u(k,h) (12) 



x u(p\\> N )^G M Y-F 2 — + F 3 ^-— j u (p,\ N ) 



Knowing both expansions of the scattering amplitude is useful, particularly when making 
comparison to other work. Our analysis will primarily use the second expansion, with the 
invariants denoted with tildes. A selection of expressions using the primed invariants is 
given in Appendix EI 

The scalar quantities Gm, F 2 , and F 3 are complex functions of two variables, say v and 
Q 2 . We will also use 

G e = G m -(1 + t)F 2 . (13) 

In order to easily identify the one- and two-photon exchange contributions, we introduce 
the notation Gm = Gm + 8Gm, and Ge = Ge + SGe, where Gm and Ge are the usual 



proton magnetic and electric form factors, which are functions of Q 2 only and are defined 
from matrix elements of the electromagnetic current. The amplitudes F3 = 5F%, 5Gm, and 
SGe, originate from processes involving the exchange of at least two photons, and are of 
order e 2 (relative to the factor e 2 in Eq. (|12jl). 
The cross section without polarization is 

da tor do tvs 



dQ Lab e(l + t) dVL Lab 

where r = Q 2 /(4M 2 ), e is 



(14) 



2 0\~ ( s - u) 2 + t(4M 2 - 1) 



l + 2(l + x)W-j = ,_„ 3 _ t4itf ,_ t ■ US) 



9 is the electron Lab scattering angle, the "no structure" cross section is 

da NS _ ^a 2 cos 2 § E' 3 



(16) 



dn Lab g 4 -E 

and E and £" are the incoming and outgoing electron Lab energies. For one-photon ex- 
change, e is the polarization parameter of the virtual photon. The reduced cross section 
including the two-photon exchange correction is given by [7] 



a R = G 2 M + £ -G 2 E + 2G M K( K 5G M + e^F^+2 £ -G E n( K 5G E + ^F^+0(e A ), (17) 

where 1Z stands for the real part. Comparison to results elsewhere is often facilitated by the 
expression 



W = Jap = V T(l + T) — ■ (l8) 

The general expressions for the double polarization observables for an electron beam of 
positive helicity (h = +1/2) and for a recoil proton polarization along its momentum (Pi) 
or perpendicular, but in the scattering plane, to its momentum (P s ) can be derived as (for 
m e = 0) Q: 

P s = A s = (19) 
= ~ \j 2£ ^~ £) ~\- {GeG m + G E TZ (8Gm) + G M 1l(6G E + ^F 3 ) + 0(e 4 )} , 

P l = -Ai = v / T 3 ^^ \g 2 m + 2G M TZ \6G M + j^^Fs) + ^(e 4 
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The polarizations are related to the analyzing powers A s or Ai by time-reversal invariance, 
as indicated above. Note that Pi is precisely unity in the backward direction, e — 0. This 
follows generally from lepton helicity conservation and angular momentum conservation. 

An observable which is directly proportional to the two- (or multi-) photon exchange is a 
single-spin observable which is given by the elastic scattering of an unpolarized electron on a 
proton target polarized normal to the scattering plane (or the recoil polarization P n normal 
to the scattering plane, which is exactly the same assuming time- reversal invariance). The 
corresponding single-spin asymmetry, which we refer to as the target (or recoil) normal spin 
asymmetry (A n ), is related to the absorptive part of the elastic eN scattering amplitude 
(27 1 . Since the one-photon exchange amplitude is purely real, the leading contribution to 
A n is of order 0(e 2 ), and is due to an interference between one- and two-photon exchange 
amplitudes. The general expression for A n in terms of the invariants for electron-nucleon 
elastic scattering is given by (in the limit m e = 0), 



P n = A n = (20) 

- T R {- G » 1 ( SG ° + W^) + Ge1 ( sg ™ + (if;) W^) } • 



T 

where X denotes the imaginary part. 

Another single-spin observable is the normal beam asymmetry, which is discussed else- 
where j3| and which is also zero in the one-photon exchange approximation. It is propor- 
tional to the electron mass, and the asymmetry is of O(10~ 6 ) for GeV electrons. It is possibly 
observable in low energy elastic muon-proton scattering. It was measured in experiments 
at MIT/Bates and MAMI at electron beam energies below 1 GeV. It is possibly also 
observable in low-energy elastic muon-proton scattering. 



III. THE TWO-PHOTON EXCHANGE CONTRIBUTION TO ELASTIC 
ELECTRON-QUARK SCATTERING 

In order to estimate the two-photon exchange contribution to Gm, F 2 and F 3 at large 
momentum transfers, we will consider a partonic calculation illustrated in Fig. ^ To begin, 
we calculate the subprocess on a quark, denoted by the scattering amplitude H in Fig. d 
Subsequently, we shall embed the quarks in the proton as described through the nucleon's 
generalized parton distributions (GPD's). 
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FIG. 1: Handbag approximation for the elastic lepton-nucleon scattering at large momentum 
transfers. In the partonic scattering process (indicated by H), the lepton scatters from quarks in 
the nucleoli, with momenta p q and p' . The lower blob represents the GPD's of the nucleon. 

Elastic lepton-quark scattering, 

l(k) + q(p q ) -> l(k') + q(p' q ) , (21) 

is described by two independent kinematical invariants, s = (k + p q ) 2 and Q 2 = —t = 
— {k — k') 2 . We also introduce the crossing variable u = (k—p' q ) 2 , which satisfies s + u = Q 2 . 
The T-matrix for the two-photon part of the electron-quark scattering can be written as 

H h , x = ^-u(k',h) llM u(k,h) ■ u(p' q ,X) (A 7 " + hl-KP^) u(p q , A), (22) 

with P q = (p q + p')/2, where e q is the fractional quark charge (for a flavor q), and where 
u(p q , A) and u(p' q , A) are the quark spinors with quark helicity A = ±1/2, which is conserved 
in the scattering process for massless quarks. Quark helicity conservation leads to the 
absence of any analog of F 2 in the general expansion of Eq. (|12j) . 

In order to calculate the partonic scattering helicity amplitudes Hh,\ of Eq. ()22|) at order 
0(e 4 ), we consider the two-photon exchange direct and crossed box diagrams of Fig. |21 
The two-photon exchange contribution to the elastic electron-scattering off spin 1/2 Dirac 
particles was first calculated in Ref. j^ij, which we verified explicitly. For further use, we 
separate the amplitude fi for the scattering of massless electrons off massless quarks into 
a soft and hard part, i.e. fi = + f^ ard . The soft part corresponds with the situation 
where one of the photons in Fig. ^carries zero four-momentum, and is obtained by replacing 
the other photon's four-momentum by q in both numerator and denominator of the loop 
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FIG. 2: Direct and crossed box diagrams to describe the two-photon exchange contribution to the 
lepton-quark scattering process, corresponding with the blob denoted by H in Fig. ^ 
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integral [31]. This yields, 
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(23) 
(24) 



where fl°^, which contains a term proportional to In A 2 (A is an infinitesimal photon mass), 
is IR divergent. The amplitude / 3 resulting from the diagrams of Fig. EJ is IR finite, and its 
real part is 

e 2 1 



K If; 



— <s\n 

Air z s u 



+ u In 



Q 2 



s — u 



4 In 2 

u 



Q 1 



-^ln 2 



it 



Q 2 



U r 

— 71 

S 



(25) 

The correction to the electron-quark elastic cross section can be obtained from Eq. (|17j). 



da = da^ 

= da Xl (1 + 5 2l ) , 



(26) 



where da^ is the cross section in the one-photon exchange approximation and e 
—2 su / (s 2 + u 2 ) in the massless limit. Using Eqs. (}2lfl I2*4ll23|) . we obtain (for e q = +1) 

,2 ( / A 2 



20 



+ 



4. 2 l 21n ^' 

(s-u)Q 2 



(27) 



2 (s 2 + u 2 



In' 



Q 1 



In 2 



Q 2 



7T 



Q 4 



s 2 + u 2 



Q 2 



In 



Q 2 



In 



Q 2 



which is in agreement with the corresponding expression for electron-muon scattering ob- 
tained in Ref. j^ . The expressions of /i and f% can also be obtained through crossing 
from the corresponding expressions of the box diagrams for the process e + e~ — > fi + n~ as 



calculated in Ref, 
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We will need the expressions for the imaginary parts of f\ and in order to calculate 
the normal spin asymmetry A n . These imaginary parts originate solely from the direct 
two-photon exchange box diagram of Fig. |21 and are given by 

i(fr") = ~>»(T). (28) 



1 (/?""') 



4tt \ s 
e 2 (Q 2 , I S \ 



i «{2u ,n {(T 2 ) + 2^ <29) 



Notice that the IR divergent part in Eq. (|28|) does not contribute when calculating the 
normal spin asymmetry A n of Eq. (f2T)|) . Indeed at the quark level, one may complete the 
calculation for quark mass m q nonzero and find that A n is given by, 



e q e 2 m q y/gT+g Q 2 (Q 2 + Am 2 q ) 
n An 2Q l+4em 2 /Q 2 s(s - u) 1 ' 



an IR finite quantity (cf. 



IV. THE HANDBAG CALCULATION OF THE TWO-PHOTON EXCHANGE 
CONTRIBUTION TO ELASTIC ELECTRON-NUCLEON SCATTERING 

Having calculated the partonic subprocess, we next discuss how to embed the quarks in 
the nucleon. We begin by discussing the soft contributions. The handbag diagrams dis- 
cussed so far have both photons coupled to the same quark. There are also contributions 
from processes where the photons interact with different quarks. One can show that the IR 
contributions from these processes, which are proportional to the products of the charges of 
the interacting quarks, added to the soft contributions from the handbag diagrams give the 
same result as the soft contributions calculated with just a nucleon intermediate state (3^- 
Thus the low energy theorem for Compton scattering is satisfied. As discussed in the in- 
troduction, the hard parts which appear when the photons couple to different quarks, the 
so-called cat's ears diagrams, are neglected in the handbag approximation. 

For the real parts, the IR divergence arising from the direct and crossed box diagrams, 
at the nucleon level, is cancelled when adding the bremsstrahlung contribution from the 
interference of diagrams where a soft photon is emitted from the electron and from the 
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proton. This provides a radiative correction term from the soft part of the boxes plus 
electron-proton bremsstrahlung which added to the lowest order term may be written as 



Vsoft - 0"1 7 (1 + $2y,soft + Cems) 



(32) 



where <Ti 7 is the one-photon exchange cross section. In Eq. 
due to the nucleon box diagram is given by 



the soft-photon contribution 



2-y, soft 




X 2 



s-M 2 )\u-M 2 
— M 2 



f n 

where L is the Spence function defined by 

L(z 




(33) 
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u - M i 
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~2 



dt 



ln(l - t) 



(34) 



, t 

The bremsstrahlung contribution where a soft photon is emitted from an electron and proton 
line (i.e., by cutting one of the (soft) photon lines in Fig.|2j) was calculated in Ref. Q], which 
we verified explicitly, and is for the case that the outgoing electron is detected, 



brems 



2ix 2 
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4 (AEf 
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X 2 y(u-M 2 f 
1 s-M 2 
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M 2 



M 2 



u 



u 



(35) 



y M 2 — u J \ y s — M 2 
E' e el — E e is the difference of the measured outgoing electron lab energy (E' e ) 



where AE 

from its elastic value (E^ el ), and y = (y/r + \/l 
of Eqs. 



T 



One indeed verifies that the sum 



is IR finite. When comparing with elastic ep cross section data, which are 
usually radiatively corrected using 



only the difference of our $2 7 , soft + ^ ems relative to the 0(Z ) part, in their notation, of the 
radiative correction in [||. Except for the 7r 2 /2 term in Eq. (J33)) . this difference was found 
to be below 10 -3 for all kinematics considered in Fig. El 

Having discussed the two-photon exchange contribution on the nucleon when one of the 
two photons is soft, we next discuss the contribution which arises from the hard part (that 
is, neither photons soft) of the partonic amplitude coming from the box diagrams. This 
part of the amplitude is calculated, in the kinematical regime where s, —u, and Q 2 are large 
compared to a hadronic scale (s, — u, Q 2 » M 2 ), as a convolution between a hard scattering 
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electron-quark amplitude and a soft matrix element on the nucleon. It is convenient to choose 
a frame where q + = 0, as in j^, where we introduce light-cone variables oc (a°±a 3 ) and 
choose the z-axis along the direction of P 3 (so that P has a large + component). We use 
the symmetric frame, as in jl^ . where the external momenta and q are 

l Q 2 l 



k 
k' 

<1 

P 
1 

V 



nP+ 1 Q2 V 



[0,0, gj , 
1 



P- 

P- 



P+ 
1 
P+ 



M 



— \M 2 + — 



4 

Q_ 2 
4 



Then, 



4?7 



Q 2 + (1 + ^)M 2 



2 g± 



77) 2 



(36) 



4?/ 



Q 2 + (1 - V )M 2 



(37) 



one may check that s + u = 2M 2 + Q 2 and also solve for the lepton light-front momentum 
fractions, 77 = k + / P + = k' + /P + , as 

1 



n 



Q 2 + AM 2 



u 



2 VM 4 - s 



(38) 



For comparison, forward scattering in the CM, 9cm = 0°, matches to 77 = and backward 
scattering, # CM = 180°, matches to r] = (s — M 2 )/s. 

In the q + = frame, the parton light-front momentum fractions are defined as x = 
Pq /P + = p' q + /P + . The active partons, on which the hard scattering takes place, are approx- 
imately on-shell. In the symmetric frame, we take the spectator partons to have transverse 
momenta that are small (relative to P) and can be neglected when evaluating the hard 
scattering amplitude H in Fig. ^ The Mandelstam variables for the process (|21j) on the 
quark, which enter in the evaluation of the hard scattering amplitude, are then given by 



(x + 77) 1 



Q 5 



u 



(x — T]f 



Q 2 



(39) 



Axr] ~* Axt] 

Note that in the limit x ~ 1, where s ~ s and u ~ u , the quark momenta are collinear with 
their parent hadron momenta, i.e. p q ~ p and p' q ~ p' . This is the simplest situation for 
the handbag approximation, in which it was shown possible to factorize the wide angle real 
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Compton scattering amplitude in terms of a hard scattering process and a soft overlap of 
hadronic light-cone wave functions, which in turn can be expressed as moments of generalized 
parton distributions (GPD 's) 0Q- In the following we will extend the handbag Q 
371 ] formalism to calculate the two-photon exchange amplitude to elastic electron-nucleon 
scattering at moderately large momentum transfers, and derive the amplitude within a 
more general unfactorized framework by keeping the x dependence in the hard scattering 
amplitude (i.e., by not taking the x — > 1 limit from the outset). 

For the process (J^J) in the kinematical regime s, —u, Q 2 » M 2 , the (unfactorized) hand- 
bag approximation implies that the T-matrix can be written as 2 

1 dx 



rphard 



i 

X - 
2 



H q (x,0,q 2 ) u{p' 1 X N ) 1 -nu{p 1 \ N ) + E« (x, 0, q 2 ) u(p', X' N ) ^ u(p, \ N ) 
- E 2 K+t - F £ -|] ■ 2 Sgm(x) H9 (X ' °' ^ X ' n) 7 ' Hl5 U{P > Ajv) ' 

1 X q 

where the hard scattering amplitude H hard is evaluated using the hard part of j\ and fz, with 
kinematics s and u according to Eq. (|39|). and where n M is a Sudakov four- vector (n 2 = 0), 
which can be expressed as 

{-77 P" + K"} . (41) 



VM 4 - su 

Furthermore in Eq. ()4U|). H q , E q , H q are the GPD's for a quark q in the nucleon (for a review 
see, e.g., Ref. [3]). 

^From Eqs. (J12)) . (|22|) . and (|4(J|) the hard 27 exchange contributions to 5Gm> 5Ge, and 
F3 are obtained (after some algebra) as 

5 Qhard = ^ ^ 

5G h E ard = - (i±£) (A - C) + ^^5, (43) 
A-^Wc). (44, 



2 The corresponding equation in Ref. 0] contains typographical errors regarding factors of (1/2). The 
remaining equations in that paper are written correctly. 
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with 



A 



B 



C 





dx 


{* 


- u)f^ ard - 






X 




(s-u) 






dx 


\s 


- u)f? ard - 






X 




(s-u) 





1 

J2e 2 q (H*-t& 



1 i7t 

-/^sgm(x) VeJ^, 
1 x , 



(45) 



where note that in Eqs. (|42|) - (j44j) . the partonic amplitude fi has its soft IR divergent part 
removed as discussed before. 

Equations (|I2J)-(|IIJ) reduce to the partonic amplitudes in the limit M — > by considering 
a quark target for which the GPD's are given by, 

H q -> 5(1 -x), 



In this limit, and using the identity 



E q -> 



# 9 -> 5(1 -x). 



sm s — m 2e 



s-u 4 1 + £ 



(46) 



(47) 



we find that 



5Fo 



hard 



M 
M 2 



(48) 



^From the integrals A, £>, and C, and the usual form factors, we can directly construct 
the observables. The cross section is 



&R — &R,soft + &R,hard 



where 



&R,hard 



;i + e) G M K (A) + y/2e(l + e)- G E K (B) + (1 - e) G M K (C) 

T 



(49) 



(50) 
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^From Eqs. (p£2j) to (|3*5)l and the discussion surrounding them, we learned that to a good 
approximation the result for the soft part can be written as 

o-R, so ft = 0\R,i7 (l + vra + 5 MT ) , (51) 

where 5 MT is the correction given in Ref. ^. Since the data is very commonly corrected 
using 0, let us define u R MT corr = 0\r/(1 + 5 A/T ). Then an accurate relationship between 
the data with Mo-Tsai corrections already included and the form factors is 

MT corr ( , £ ri2 



a Ml corr = ^ + (1 + ^ + ? (53) 

where the extra terms on the right-hand-side come from two-photon exchange and 0(e 4 ) 
terms are not included. The reader may marginally improve the expression by including with 
the (1 + ira) factor the circa 0.1% difference between our actual soft results and those of jsj]; 
from our side the relevant formulas are the aforementioned to (|35|). Since the Mo-Tsai 
corrections are so commonly made in experimental papers before reporting the data, the 
"MT corr" superscript will be understood rather than explicit when we show cross section 
plots below. Finally, before discussing polarization, the fact that a ir 2 /2 term, or (na) term 



after multiplying in the overall 
criterion we used, that of Ref. 



actors, sits in the soft corrections has to do with the specific 
~l 

311 ] . to separate the soft from hard parts. The term cannot 
be eliminated; with a different criterion, however, that term can move into the hard part. 
The double polarization observables of Eqs. ()53I54|) are given by 



2g(1 r £) ^ I G E G M + G E n (C) + G M Ti (B) + 0(e A ) } , (53) 



P l = VT^— {G 2 M + G M 1l(A + C)+0(e A )} J (54) 
and the target normal spin asymmetry of Eq. (|20[) is 



A n = ]/ 2£{1 t +£) ^ \g e X(A) - yjl±±G M I{B) \ , (55) 

One sees from Eq. (|55|) that A n does not depend on the GPD H. 

We will need to specify a model for the GPD's in order to estimate the crucial integrals 
Eqs. (I45|) for the two-photon exchange amplitudes We will present results from two different 
GPD models: a gaussian model and a modified Regge model. 
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First, following Ref. |20j , we use a gaussian valence model which is unfactorized in x and 
Q 2 for the GPD's H and H, 

H*(x, 0, q 2 ) = q v (x) exp ( J±Z^- ) , (56) 

H«(x, 0, q 2 ) = Aq v (x) exp ( J 1 '^^ ) , (57) 

where q v (x) is the valence quark distribution and Aq v (x) the polarized valence quark dis- 
tribution. In the following estimates we take the unpolarized parton distributions at input 
scale Ql = 1 GeV 2 from the MRST2002 global NNLO fit Q as 

u v = 0.262 x" a69 (l - x) 3 - 50 (1 + 3.83 x a5 + 37.65 x) , 
d v = 0.061 a;- a65 (l - x) 4 - 03 (1 + 49.05 a; ' 5 + 8.65 a;) . 

For the polarized parton distributions, we adopt the recent NLO analysis of Ref. j^j], which 
at input scale Qq = 1 GeV 2 yields 

Au v = 0.505x~ a33 (l-x) 3 - 428 (l + 2.179x a5 + 14.57x), 
Ad v = -0.0185a+ a73 (l - x) 3 ' 864 (l + 35.47x a5 + 28.97x) . 

For the GPD E, whose forward limit is unknown, we adopt a valence parametrization 



multiplied with (1 — x) 2 to be consistent with the x — > 1 limit 4l|. This gives 



E\x^q 2 ) = ^-(1 -x)V(x)exp ( - (1 - ^ \ , (58) 



Ni y ' " v ' r \ Axa 
where the normalization factors N u = 1.377 and N d = 0.7554 are chosen in such a way 
that the first moments of E u and E d at Q 2 = yield the anomalous magnetic moments 
k u = 2n p + n n = 1.673 and K d = K p + 2n n = —2.033 respectively. Furthermore, the parameter 
a in Eqs. (|56I57I58|) is related to the average transverse momentum of the quarks inside the 
nucleon by a = 5 < k\ >. Its value has been estimated in Ref. |3| as a ~ 0.8 GeV 2 , which 
we will adopt in the following calculations. 

The GPD's just described were used in our shorter note jl^ . Recently, GPD's whose 
first moments give a better account of the nucleon form factors have become available 421. 

n 

These GPD's we refer to as a modified Regge model |42J, and entail 

H q (x,0,q 2 ) = q v (x) x < {1 ~ x)q2 , 
E\x^q 2 ) = ^-{l-xf"q v {x)x a '^^\ 

H q (x,0,q 2 ) = Aq v (x) x d '^- x)Q2 . (59) 
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We still use the same q v and the same Aq v as given above. The five parameters are 

a[ = 1.098 GeV" 2 , a' 2 = 1.158 GeV" 2 , a[ = 1.000 GeV~ 2 , 

r) u = 1.52, rj d = 0.31, (60) 

and the normalization factors here become N u = 1.519 and Nj, = 0.9447. The modified 
Regge GPD's formally do not give convergence at low Q 2 for integrands with negative 
powers of x, such as we have here (or as one finds injl9|). The integrals could be defined 
by analytically continuing in the Regge intercept 3> Q]. We will use them only for Q 2 > 2 
GeV 2 , and all the integrals converge straightforwardly. 

We shall investigate in forthcoming plots the sensitivity of the results to the two GPD's. 

V. RESULTS 

A. Cross section 

In Fig. El we display the effect of 2j exchange on the reduced cross section or, as given 
in Eq. (|17|). for electron-proton scattering. For the form factor ratio, we always use G P E /G P M 
as extracted from the polarization transfer experiments 3| • 

We should remind the reader that G P M is also obtained from the reduced cross section 
data: the normalization gives G P M and the slope gives G P E /G P M . As a starting point we adopt 
the parametrization for G P M , of Ref. The straight dotted curves of Fig. El show that 

the values of G P E /G P M extracted from the polarization experiments are inconsistent with the 
one-photon exchange analysis of the Rosenbluth data, corrected with just the classic Mo 
and Tsai radiative corrections jsj], in the Q 2 range where data from both methods exist. We 
then include the 27 exchange correction, using the GPD based calculation described in this 



3 



and 



paper. The plots show the results from both the GPD's used in our shorter note 
recorded in the previous section, as well as from the alternative GPD's also described in the 
previous section. The results are rather similar. 

It is also important to note the non-linearity in the Rosenbluth plot, particularly at 
the largest e values. One sees that over most of the e range, the overall slope has become 
steeper, in agreement with the experimental data. This change in slope is crucial: we see that 
including the 27 exchange allows one to reconcile the polarization transfer and Rosenbluth 
data. 
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FIG. 3: Rosenbluth plots for elastic ep scattering: cjr divided by (^l p Gd) 2 , with Gd 



J 1 + 

nn 



(5 2 /0.71)~ 2 . Dotted curves: Born approximation using Gep/Gmp from polarization data [lj, y]. 
Solid curves: full calculation using the modified Regge GPD, for the kinematical range — u > M 2 . 
Dashed curves: same as solid curves but using the gaussian GPD. The data are from Ref. |fj. 

It is clearly worthwhile to do a global re-analysis of all large Q 2 elastic data including 
the 27 exchange correction in order to redetermine the values of G P E and G P M . For example, 
in order to best fit the data when including the 27 exchange correction, one should slightly 
change the value of G P M of Ref. 3|. A full analysis is beyond the scope of this paper, 

In Fig. HI we show a similar plot for electron-neutron elastic scattering. Because of 
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Cross section for en elastic scattering 
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FIG. 4: A Rosenbluth plot for elastic en scattering. The curves with the two-photon corrections 
are plotted only for — u > M 2 . For and G\ we used fits from 



44( and 



451 ] . respectively. There 



is little slope change for the neutron case, for reasons noted in the text. 



a partial cancellation between contributions proportional to G% and G 7 ^, there is little e 
dependence in the corrections, and the slope is not appreciably modified. We took G 1 }, from 

n n 

the fit of [44J; for G% we used the fit given in [45J. 



B. Single spin asymmetry 

The single spin asymmetry A n or P n is a direct measure of the imaginary part of the 27 
exchange amplitudes. Our GPD estimate for A n for the proton is shown in the left-hand 
plot of Fig. as a function of the CM scattering angle for fixed incoming electron lab energy, 
taken here as 6 GeV. Also shown is a calculation of A n including the elastic intermediate 
state only j^. The result, which is nearly the same for either of the two GPD's that we 
use, is of order 1%. 

Fig. El on the right also shows a similar plot of the single spin asymmetry for a neutron 
target. The predicted asymmetry is of opposite sign, reflecting that the numerically largest 
term is the one proportional to Gm- The results are again of order 1% in magnitude, though 
somewhat larger for the neutron than for the proton. 

A precision measurement of A n is planned at JLab on a polarized 3 He target; it 
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FIG. 5: Nucleon analyzing power, which is equal to normal recoil polarization. The elastic con- 
tribution (nucleon intermediate state in the two-photon exchange box diagram) is shown by the 
dotted curve j^]. The GPD calculation for the inelastic contribution is shown by the dashed curve 
for the gaussian GPD, and by the solid curve for the modified Regge GPD. The GPD calculation is 
cut off in the backward direction at —u = M 2 . In the forward direction the modified Regge GPD 
result goes down to Q 2 = 2 GeV 2 and the gaussian GPD result to Q 2 = M 2 . 

will provide access to the elastic electron-neutron single-spin asymmetry from two-photon 
exchange. 

C. Polarization transfers 

The polarization transfer method for measuring the ratio Ge/Gm depends on measuring 
outgoing nucleon polarizations Pi and P s for polarized incoming electrons. Their ratio is 



Ps 

Pi 



2e G, 



r(l + e) G 



Al 



(61) 



in the one-photon exchange calculation. This also is subject to additional corrections from 
two-photon exchange. However, the impact of the corrections upon Ge is not in any way 
enhanced, and so one expects and finds that the corrections to Gg measured this way are 
smaller than the corrections to coming from the cross section experiments. 

Figure El shows in the upper two panels the calculated Pi and P s for ep scattering with 
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FIG. 6: Recoil proton polarization components P s and Pi and their ratios relative to the I7 
exchange results (lower panels) for elastic ep scattering at Q 2 = 5 GeV 2 . The dotted curves in the 
upper panels are the Born approximation (I7 exchange) results. The solid curves include the 27 
exchange correction using the GPD calculation, for the kinematical range where both s, —u > M 2 . 



and without the two-photon exchange terms, for 100% right-handed electron polarization 
and with fixed momentum transfer Q 2 = 5 GeV 2 . The two GPD's were presented in the 
previous section, and we use again the polarization Ge/Gm from "\ and Gm from J^. The 
corrections to the longitudinal polarization are quite small, as is seen again in the lower left 
panel, where the ratio of the full calculation divided by the one-photon exchange calculation 
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FIG. 7: Recoil neutron polarization components P s and Pi and their ratios relative to the I7 
exchange results (lower panels) for elastic ep scattering at Q 2 = 4.3 GeV 2 . See Fig. El for notation. 



is shown. The lower right panel shows the corrections to the P s /Pi ratio, given as a ratio 
again of the full calculation to the one-photon calculation. An experiment to measure the 
e-dependence of P s /Pi is planned at JLab j^. This will allow a test of the two-photon 
corrections. 

Figure [3 shows the corresponding plots for the neutron, at a momentum transfer squared 
of 4.3 GeV 2 . If one needs to choose between the GPD's, the modified Regge model should 
be chosen as it gives the better account of the existing data on the form factors, the neutron 
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form factors in particular 
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D. Positron-proton vs. electron-proton 

Positron-proton and electron-proton scattering have the opposite sign for the two-photon 
corrections relative to the one-photon terms. Hence one expects e + p and e~p elastic scat- 
tering to differ by a few percent. Figure |H1 shows our results for three different Q 2 values. 
These curves are obtained by adding our two-photon box calculation, minus the correspond- 
ing part of the soft only calculation in jjj, to the one-photon calculations; hence, they are 
meant to be compared to data where the corrections given in [^J have already been made. 
Each curve is based on the gaussian GPD and is cut off at low e when — u = M 2 . Early data 
from SLAC are available 48[; more precise data are anticipated from JLab Q (Ref. Q 
used the Meister-Yennie 9( soft corrections rather than those of Mo and Tsai. We have 
checked that for these kinematics the difference between them is smaller than 0.1%, which 
is negligible compared to the size of the error bars.) 



E. Possibilities at lower \u\ 

In numerical calculations, we used a conservative requirement that the values of the 
Mandelstam variable \u\ > M 2 in order to apply the partonic description. In a 'handbag' 
mechanism of wide-angle Compton scattering on a proton, such a requirement is needed 
to enforce high virtuality of the quark line between the two currents, making sure that 
short light-cone distances dominate. However, our case of electron-quark scattering via 
two-photon exchange involves 4-dimensional loop integration, and small values of \u\ do not 
necessarily mean that the struck quark has small virtuality. Analyzing the two-photon- 
exchange loop integral in terms of Sudakov variables one may show that for the backward 
(u — ► 0) electron-quark scattering, high virtuality of the quark dominates the loop integral, 
thereby justifying extension of our approach to the region of small u, as long as s and 
—t remain large. Such an analysis may be found in the literature for the backward-angle 
electron- muon scattering in QED jscj, and we found our formalism consistent with these 
early calculations. 
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FIG. 8: Ratio of e + /e~ elastic cross sections on the proton. The GPD calculations for the 27 
exchange correction are for three fixed Q 2 values of 2, 5, and 9 GeV 2 , for the kinematical range 
where — u is above M 2 . Also shown are all known data, from 48 1, with Q 2 above 1.5 GeV 2 (the 
missing central value is at 1.111). The numbers near the data give Q 2 for that point in GeV 2 . 

F. Rosenbluth determinations of Ge/Gai including 2-photon corrections 

Previous Rosenbluth determinations of Ge/Gm were made using data which had been 
radiatively corrected using the Mo-Tsai [^] or comparable [l| prescription. Given the work 
in this paper, we would now say that these corrections are just a part of the total radiative 
correction. One should also include the hard two-photon corrections. 

We present here new Rosenbluth determinations of G P E /G P AI using known data but in- 
cluding the two-photon corrections. We used cross section data from Andivahis et al. j^|, 
and made a \ 2 fit to the data at each of the five Q 2 selected using our full calculation and 
allowing both G P M and G P E /G P M to vary. We included the lowest e points in the data by 
making a linear extrapolation of our calculations from higher e. (For the record, and for the 
e's in question and to the precision we need, the result is numerically the same as doing our 
GPD calculation at these e's, even though \u\ is below M 2 .) 
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FIG. 9: Rosenbluth determinations of Ge/Gm including the 2-photon corrections. The polarization 



data is from Jones et al. and Gayou et al. [2|], and the Rosenbluth determinations without the 
two-photon corrections are from 5l(. The Rosenbluth Ge/Gm are based on data from Andivahis 



et al. [g]. Some of our points for the Rosenbluth results are slightly offset horizontally for clarity. 

The results are shown in Fig. El The figure also shows the results of the polarization 
transfer measurements, and Rosenbluth results taken from j^jj, which do not include the 
hard two-photon corrections. The polarization results also have radiative corrections, but 
the size of them is, as one has learned from Fig. El smaller than the dots of the data points. 
The solid squares in Fig. El show the Ge/Gm ratios we have extracted with Ref. [f| data and 
the two-photon corrections with the gaussian GPD. The results with the modified Regge 
GPD are omitted to reduce clutter on the graph; they are about the same as for the gaussian 
for Q 2 of 2-3 GeV 2 , and a bit larger at the higher Q 2 . 

For Q 2 in the 2-3 GeV 2 range, the Ge/Gm extracted using the Rosenbluth method 
including the two-photon corrections agree well with the polarization transfer results. At 
higher Q 2 , there is at least partial reconciliation between the two methods. 

One may comment on the growth of the error bars at higher Q 2 . The calculation with 
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the two-photon contributions includes a lowest order term quadratic in G P E and a correction 
linear in G P E with opposite sign. The partial cancellation explains the reduced sensitivity to 
changes in G E . 

VI. CONCLUSIONS 

We have studied the effects of two-photon physics for lepton-nucleon elastic scattering. 
Our main result is a calculation of the two-photon exchange contributions including contri- 
butions coming when intermediate particles which are far off shell. The main impediment 
to performing this calculation is the lack of knowledge of nucleon structure. Here we have 
used a partonic "handbag" model to express the contributions when both photons are hard 
in terms of the generalized parton distributions (GPD's) of the nucleon. The GPD's also en- 
ter calculations of deeply virtual Compton scattering, wide angle Compton scattering, and 
exclusive meson photoproduction, which are consistent with models for the GPD's. The 
calculations which we have presented are valid when s, —u, and Q 2 are large, although we 
have argued in subsection IV El and Ref . , that the requirement on — u is not compulsory 
for eiV elastic scattering). We have presented our results requiring that the magnitude of 
each of the invariants is above M 2 . 

We have found that in Rosenbluth plots of the differential cross section vs. e, that the 
two-photon exchange corrections gives an additional slope which is sufficient to reconcile 
qualitatively the difference between the Rosenbluth and polarization data. The change in 
the effective slope in the Rosenbluth plots comes only from corrections where both photons 
are hard. The reconciliation thus implies only a minor change in the Ge/Gm ratio as 
obtained from the polarization data, since those data receive smaller two-photon corrections 
to Ge/Gm- 

Two-photon exchange has additional consequences which could be experimentally ob- 
served. For polarizations P s and Pi, there are two-photon corrections which are small but 
measurable. For the normal direction, the polarization or analyzing power is zero in the one- 
photon exchange limit, but the presence of the two-photon exchange amplitude leads to a 
nonzero effect of (9(1%). We also predict a 0(iew%) positron-proton/electron-proton asym- 
metry. The predicted Rosenbluth plot is no longer precisely linear; it acquires a measurable 
curvature, particularly at high e. 
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Thus, in summary, we have shown that the hard two-photon exchange mechanism sub- 
stantially reconciles the Rosenbluth and polarization transfer measurements of the proton 
electromagnetic elastic form factors. We have also emphasized that there are important 
experimentally testable consequences of the two-photon amplitude. 
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APPENDIX A: CROSS SECTION AND POLARIZATION RESULTS IN THE 
AXIAL- VECTOR REPRESENTATION 

This appendix records the cross section and polarization results using the expansion of 
the scattering amplitude in the axial- vector representation given by Eq. 



,2 



T h ,\' N \ N = { u(k', h)^^u(k, h) x u(p', X' N ) 



7 "G" M F' 

1 M M 2 



u(p,X 



+ u(k , h)-f^ 5 u(k, h) x u(p, X N ) ffG^ u(p, Ajy) | • (Al) 

1. Form factors and observables 

The scalar invariants or form factors are in general complex and functions of two variables. 
We also define 

G' E = G' M - (1 + r)F 2 ' . (A2) 

The relations between the present scalar invariants and the ones used in most of the text 
follow from Eq. (jllj) and are 

G' M = G M + ^T^3 



F' 2 = F 2 
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G' A = -tF 3 

G'e = G E + ^F 3 . (A3) 

The invariants may be separated into parts coming from one-photon exchange and parts 
from two- or more-photon exchange, 

G'm = Gm + SG' M , 
G'e — Ge + SG' E , 

G'a = *G' A , (A4) 

where Gm(Q 2 ) and Ge{Q 2 ) are the usual magnetic and electric form factors, defined from 
matrix elements of the electromagnetic current and real for spacelike Q 2 . The quantities 
SG' M , SG' E , and G' A are C(e 2 ) relative to Gm or Ge- 
The reduced cross section in Eq. ((T^) is 



°R = \G' M \ 2 + £ -\G' E \ 2 + (1 + T y £2) G M K (G' a ) + 0(e A ) . (A5) 
The polarizations of the outgoing nucleons or analyzing powers of the target nucleons are 



A n = \l v r ' ' — I T{G%G' M ) + ^-—.-—G E 1(G' a ) + 0(e 



2e(l + e) 1 x , /1 + r 1-e ^, 4 . 



A s = - P e \l v r ; — { K (G%G' M ) + x l—.-—G E n (G' A ) + 0(e 



2e(l-e) 1 I _ / ^,*^ / s , /1 + r 1 + £ ~ , ~, ^ 



P = -A = P e ^|v / l^^|G';, / | 2 + 2^^G, / ^(G^) + 0(e 4 )| . (A6) 

The only single spin asymmetry is P n or A n . Further, P n or A n is zero if there be only 
one-photon exchange, so observation of a non-zero value is definitive evidence for multiple- 
photon exchange. Polarizations P s or Pi are double polarizations. The expressions for them 
are proportional to the electron longitudinal polarization P e (with, e.g., P e —li£h — +1/2). 



2. Electron-quark elastic scattering amplitudes 

The two-photon part of electron-quark elastic scattering is given by 

lee) 2 ( _ , 
H hj x = I < g M u(k', fe)7 M tt(fe, h) ■ u(p' q , X)j^u(p q , A) 
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(27) - 
9 A u 



ee„ 



Q 2 



(k',h)^ 5 u(k,h) ■ u(p' q , \)Yl 5 u(p q ,\)} 
s-u- {2h)(2\)tj g M + ((2/i)(2A)(a - u) - tjg™} , 



(A7) 



when the electrons and quarks are both massless. The theorem of Eq. (|11|) relates 



9m — Ji H : — J3 > 



(27) 
9 a 



4/3 • 



(A* 



We split the two-photon part of into a hard and soft part, = g^ 1 + g^ ard , using 



e jsijl, and have 
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3. Embedding 



The nucleon form factors are given in terms of the quark-level amplitudes and generalized 
parton distributions by 



r^/(hard) 
0Lr M 



l + £ -A- l -^C 



2s 



2s 



5G 



/(hard) 



5G 



/(hard) 



2e 

t l + e 



(A-C) 



s — u 2s 

Quantities A, B, and C are the same as in the text, but now written as 

<2 



.4 



X 



s — u 



(A10) 
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J X S 11 

sgm(x) ^ tf 9 , 



(All) 



and it is understood in Eq. (|A10|) that the partonic amplitude #m has its soft part removed. 



APPENDIX B: QUARK MASS SENSITIVITY 

1. Kinematics, and imaginary parts of the hard amplitudes 

We have until now set the quark mass to zero. 

To investigate how severe this approximation is, we will examine the effect of restoring 
the quark mass for the analyzing power calculations, though still keeping only the quark 
chirality conserving amplitudes. There are three modifications. The expressions for s and u 
become 



(x + rfy 



Q 2 + 



X + X] 



(X 



-m~ 



u 



Q 2 + 



x-r) 2 



x 



(Bl) 



Axt] x q AXT] 

where m q is the effective quark mass. The general electron-quark scattering amplitude, 
Eq. (|22|) ) should have another term with a scalar function that we may call J2 in analogy 
with the expansion of the electron-nucleon amplitude given in Eq. ()12j1 . However, this term 
flips quark helicities, and presently the formalism for embedding quark amplitudes into the 
nucleon using GPD's involves only the non-chirality flip GPD's. There is neither theoretical 
development nor experimental information regarding chirality flip GPD's, and so we shall 
ignore f% as well as helicity flip parts of other amplitudes. Including the quark mass leads 
to a modification of the hard scattering amplitudes so that 



Tjhard , Tjhard 



= q5 \i S ~ U ] fl ~ SU f3 



Q~ 2 



2m 2 

u-^- 2 Q 2 

s — mi 



fi hard + IK 



su 
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(B2) 



within the quantities A and B (C is not needed for the analyzing power), and 
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FIG. 10: Quark mass correction plots for the analyzing power or normal polarization, for fixed 
electron incoming lab energies of 3 and 6 GeV. 
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sQ 2 — (s — m 2 y 

The results for the analyzing power A n when including quark masses in the quark helicity 
conserving amplitudes are shown in Fig. f° r quark masses 300 MeV and 450 MeV. The 
effects are clearly not large. 



2. Real parts of hard two-photon exchange amplitudes with finite quark mass 

When the quark mass is not zero, we have 
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where s and u were given in Eq. (IB II) . and 
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The effect of the quark mass corrections upon the reduced cross section is shown in Fig. ^2 
for Q 2 = 4 GeV 2 . One sees from Fig. ^2 that the quark mass effects mainly influence our 
result at small values of e, where \u\ becomes small. They show the theoretical error on our 
calculation in this region. A full calculation also requires quantifying the effect of the cat's 
ears diagrams. A study of such corrections is clearly worthwhile for a future work, both for 
two-photon exchange amplitudes and for wide-angle Compton scattering. 
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